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THEORETICAL  AND  EXPERIMENTAL  STUDIES 
RELATING  TO  THE  PURGING  OF  RADIOACTIVITY 
FROM  A  GAS  WELL  STIMULATED 
BY  A  NUCLEAR  EXPLOSION 

*  Introduction 

One  question  associated  with  the  use  of  nuclear  explosions  for  natural  gas  well 

stimulation  is  that  of  eliminating  or  controlling  the  radioactivity  in  the  gas  produced 

❖ 

from  the  well.  A  proposed  solution  is  to  purge  the  chimney  before  permitting  the  gas 
to  enter  the  pipeline,  and  to  phase-in  the  production  of  the  well  so  as  to  obtain  a  dilu¬ 
tion  of  activity  to  an  acceptable  level.  The  variation  in  the  radioactivity  concentration 
during  a  purge  is  pertinent  to  the  problem.  A  previous  analysis^  was  reviewed  and  it 
was  determined  that  further  study  on  this  problem  was  necessary.  After  being  revised 
and  extended,  the  resultant  estimated  variation  of  radioactivity  concentrations  during 
well  production  for  a  range  of  chimney  geometries  was  prepared  and  compared  with 
uniform  mixing  models.  The  revised  analysis  indicates  that  the  quantity  of  gas  which 
may  be  flared  to  remove  substantially  all  the  radioactivities  may  be  only  about  one-half 
that  previously  estimated. 

*  The  classic  uniform  mixing  model  used  to  evaluate  mixing  in  a  large  tank  of 
constant  in-and-out  flow  leads  to  continuous  exponential  decrease  in  radioactivity  con- 

^  centration  which  approaches  but  never  reaches  zero.  The  displacement  model  pre¬ 

sented  here  leads  to  a  small  volume  of  release  at  constant  concentration  followed  by  a 
rapid  decrease,  finally  going  to  zero  in  two  to  four  chimney  volumes.  The  concentra¬ 
tion  deduced  from  the  displacement  model  is  initially  higher  but  at  late  times  less  than 
that  deduced  from  the  constant  mixing  model.  Experimental  tests  of  both  models  are 
discussed. 

Theoretical  Analysis 

First,  assume  a  nuclear  chimney  in  the  idealized  form  of  a  vertical,  right 
circular  cylinder  with  radius,  R,  and  height,  H.  Also  assume  that  this  chimney  is 
located  in  a  natural  gas  reservoir,  and  that  the  top  and  bottom  chimney  surfaces 
coincide  with  the  top  and  bottom  of  the  reservoir.  The  reservoir  may  be  of  infinite 


For  a  description  of  the  nuclear  ''chimney*'  and  other  features  produced  by  under¬ 
ground  nuclear  explosions,  see  UCRL-14756,  "Industrial  Applications  of  Contained 
Nuclear  Explosions,"  D.  E.  Rawson,  July  1966, 
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horizontal  extent,  or  it  xnay  be  finite,,  in  the  form  of  a  right  circular  cylinder  with  a 
radius  much  greater  than  the  chimney  radius,  R.  The  top  and  bottom  surfaces  of  the 
chimney  and  reservoir  are  impermeable  except  for  a  non -penetrating  drainage  well 
located  at  the  top  center  of  the  chimney.  The  chimney  properties  such  as  porosity 
and  permeability  are  uniform  and  isotropic,  as  are  the  properties  of  the  reservoir 
outside  the  chimney.  The  initial  values  of  gas  pressure  and  temperature  are  uniform 
throughout  the  chimney  and  reservoir,  and  the  gas  itself  is  an  ideal  gas  of  constant 
composition.  Therefore,  a  radially  symmetric  system  exists  in  which  the  distances, 
r,  from  the  chimney  centerline  and  z,  from  the  top  of  the  chimney,  are  the  only  spatial 
coordinates  that  need  be  considered. 

Now,  consider  the  case  in  which  the  chimney  permeability  is  on  the  order  of 
megadarcys  and  the  reservoir  permeability  is  on  the  order  of  millidarcys.  The  model 
implicitly  assumes  the  fractured  zone  with  intermediate  permeabilities  has  zero 
thickness.  Then,  decrease  the  pressure  at  the  top  center  drainage  well  suddenly  by 
some  small  amount.  The  resultant  rarefaction  wave  will  be  rapidly  propagated 
throughout  the  chimney,  but  its  velocity  beyond  the  chimney-reservoir  interface  will 
be  relatively  imperceptible.  Similarly,  the  vertical  pressure  gradient  inside  the 
chimney  along  the  chimney-reservoir  interface  (0p/8z)  will  be  negligible  compared  to 
the  horizontal  gradient  in  the  reservoir  just  outside  the  interface  (0p/0r),  This 
permits  the  approximation  that  the  perturbed  pressure  is  uniform  over  the  entire 
interface  (if  gravity  is  neglected).  Also  assume  that  the  pressure  at  the  drainage  well 
remains  at  a  value  slightly  less  than  that  at  the  interface.  Therefore  gas  flows  from 
the  reservoir  into  the  chimney  and  up  the  well  at  a  constant  rate.  Heat  transfer  from 
the  porous  material  in  the  chimney  and  reservoir  maintains  isothermal  flow.  The  gas 
flow  in  the  chimney  is  essentially  incompressible  if  the  chimney  pressure  differential 
is  sufficiently  small. 

All  assumptions  are  compatible  with  the  existence  of  Darcy  flow  within  the 
reservoir  and  chimney.  The  flow  velocity,  v,  at  a  point  along  a  streamline  is  there¬ 
fore  given  by  the  relation  (again  neglecting  gravity): 


Ji  ^ 

0S 


(1) 


where  K  is  the  permeability,  /u  the  viscosity,  cf)  the  porosity  and  (0p/ 0s)  the  local 
pressure  gradient  along  the  streamline. 

Next,  assume  that  the  reservoir  pressure  is  uniform,  and  that  the  pressure  may 
be  considered  uniform  over  the  entire  reservoir-chimney  interface.  Then  the  pressure 
distribution  in  the  reservoir  outside  the  chimney  is  independent  of  the  coordinate,  z, 
and  is  a  function  of  only  r  and  t.  Consequently,  the  uniform  properties  of  the  reservoir, 
and  Eq,  (1),  v  is  independent  of  z  and  is  a  function  of  only  r  and  t  within  the  reservoir. 
Therefore  the  gas  velocity  entering  the  chimney,  v^,  is  uniform  over  the  entire 


-2- 


chimney-reservoir  interface.  For  isothermal,  approximately  incompressible  flow, 
this  means  that  the  flux  density  entering  the  chimney  is  similarly  uniform.  Therefore 
the  volumetric  flow  rate  into  the  chimney  through  an  elemental  area  with  height, 
dz  is; 

* 

dq^  =  (25rRdz)v^  ,  (2) 

through  the  interval  (0  <  z  <  z^)  is 

q.  =  27rRZiV.  ,  (3) 

u  1  1 

and  into  the  entire  chimney  is 


Q.  =  27rRHv.  .  (4) 

For  steady,  incompressible  flow,  Eq.  (4)  also  gives  the  solution  for  the  total 
flow  rate  out  of  the  chimney.  The  direct  proportional  relation  between  and  H  might 
appear  to  differ  from  the  result  for  the  production  rate  of  partially  penetrating  wells 
in  uniform  reservoirs  (which  corresponds  to  the  uniform  chimney).  Beyond  the  first 
few  feet  of  sand  (below  the  well  bottom)  the  additional  layers  of  sand  give  successively 
decreasing  contributions  to  the  production  from  the  well.  This  effect  is  associated 

*  with  the  nonuniform  flux  density  at  the  well  inlet  (see  Ref.  2,  p.  269-270  and  27  3). 
Furthermore,  the  solution  (Ref.  2,  p.  270-271  and  Ref.  3)  shows  a  transition  from 
spherical  flow  (with  nonuniform  flux  distribution)  in  the  vicinity  of  a  "nonpenetrating" 
well  (a  chimney)  to  radial  flow  (with  uniform  flux  density)  near  the  external  boundary 
(chimney  reservoir  interface).  Since  H  is  generally  greater  than  R,  the  constant 
pressure  interface  imposes  the  uniform  radial  flow  condition  relatively  close  to  the 
well. 

If  the  solution  for  a  short  disk  is  applied  to  a  tall  cylinder,  the  length  of  the 
streamlines  within  the  chimney  is  such  that  L(0)  =R,  L(H)  =R  +  HandL(z)  R  +  z.  Thus 
the  approximation: 

L(z)  =  R  +  z.  (5) 

^  can  be  made. 

Assuming  that  the  transit  time  for  a  rarefaction  wave  to  travel  throughout  the 

•  chimney  is  approximately  zero  (it  is  exactly  zero  for  incompressible  flow),  and  that 
the  pressure  differential  Ap  between  the  chimney  drainage  well  and  the  chimney- 
reservoir  interface  is  kept  at  a  constant  value,  then  Eq.  (1)  may  be  integrated  and 
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combined  with  Eq,  (5)  to  give  the  following  relation  for  the  average  velocity  v(z)  along 
a  given  streamline  that  enters  the  chimney  at  z: 


v(z) 


Yc  L<z)  1 

\  (9p/9s)ds 

i  \ 

-'o 

=  -  JL 

/  \ 

L(z) 

\R  +  zj 

(6) 


It  must  be  noted  that  this  is  the  spatial  average  velocity  which  is  not,  in  general,  equal 
to  the  usual  definition  of  average  velocity:  Distance  divided  by  transit  time.  This 
assumption  is  more  crude  than  that  implicit  in  Eq.  (5).  Nevertheless,  in  the  interest 
of  avoiding  complex  mathematics  like  that  in  Refs.  2  and  3,  the  transit  time  along  a 
given  streamline  through  the  chimney  can  be  defined  as: 


t(z)  =  L(z)/v(z). 


(7) 


Then  the  transit  time  along  the  z  =  0  streamline  is  a  xinit  reference  time.  It  follows 
that  Eqs.  (5)  through  (7)  may  be  combined  to  give  a  dimensionless  time: 


At  zero  time  the  chimney  is  filled  with  radioactive  gas  and  the  gas  in  the 
reservoir  is  free  of  radioactivity.  If  there  is  no  diffusion  or  dispersion  during  the 
purging  process,  purging  is  a  simple  function  of  fluid  displacement.  In  general,  as 
radioactive  gas  is  being  purged  from  the  chimney,  radioactive  gas  will  flow  through 
one  portion  of  the  drainage  well  area  Aq  and  "clean  gas"  will  flow  through  the  re¬ 
maining  area  A^: 

A^(t)  +  A^(t)  =  Aq.  (9) 


Let  z,  be  the  coordinate  of  the  streamline  that  forms  the  boundary  between  the  purged 
and  unpurged  portions  of  the  chimney.  Then  from  Eqs.  (8)  and  (9), 


and 


for 

0  <  t*  <  1, 

(10a) 

<  H  for 

1  <  t*  <  [  1  +  (H/R) 

aob) 

for 

t*  =  [1  +  (H/R)]^. 

(10c) 

If  A^/Aq  is  assumed  equal  to  z^/H,  the  relation  used  by  Higgins  and  Rodean  to 
plot  Fig.  7  of  Ref.  1  is  A^/Aq  =  1  -  z^/H,  However,  A^/Aq  f  unless  the  well 

flux  distribution  is  uniform  (which  it  is  not  as  noted  above). 

A  different  approach  that  recognizes  nonuniform  flux  distribution  at  the  drainage 
well  utilizes  Eq.  (6).  The  average  of  v(z)  for  all  the  streamlines  becomes: 
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^(H)  .  i 


'H  _ 
v(z)dz 

0 


=  KAp 
“  “m4>h 


(11) 


Similarly,  the  average  of  v(z)  for  the  unpurged  portion  of  the  chimney  is: 


v(Zj) 


fH  . 

\  v(z)(iz  =  - 


KAp 


-  z^) 


(12) 


Assuming  that  v(H)  and  viz^)  are  directly  proportional  by  the  same  factor  a  to  the 
average  velocities  through  and  A^,  respectively,  the  steady  incompressible  flow 
ratio  are: 


Q.  =  aA^  v(H) 

Qi  -  qjL  =  a  v(zp. 

Then  Eqs.  (3),  (4)  and  (11)  through  (14)  may  be  combined  to  give: 

Ir,  /  R+  H  \ 


(13) 

(14) 


r  ^ 


In 


W)  ■ 


(15) 


If  the  concentration  of  radioactivity  in  the  gas  flowing  through  is  defined  as 
imity,  then  the  mean  concentration  over  the  total  well  area  Aq  is,  from  Eq.  (15), 

/R  +  H  \ 


A 

C(z,)  =T^  = 


In 


-PF)  ■ 


(16) 


This  gives  C(z^)  =  1  when  =  0  and  C(z^)  =  0  when  z^  =  H,  Now,  from  Eqs.  (10a) 
through  (10c)  and  (16),  if  the  concentration  of  radioactivity  in  the  gas  leaving  the 
chimney  is  expressed  as  a  function  of  time. 


C(t*)  =  1 


C(t*) 


in/'-^ 

\R^/F  ) 


in 


/R  +  H\ 
[  R  j 


for 

0  <  t^'  < 

1 

(17a) 

for 

1  <  t*  < 

(-If- 

(17b) 

Then  from  Eqs.  (8),  (17a)  and  (17b),  the  average  radioactive  concentration  during  the 
entire  purging  process  can  be  found: 


(Ai)’  S. 


--{••si 


C(t*)dt*  = 


=  f-— 1 

\R  +  hJ 


(18) 
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To  make  the  integration  easier,  introduce  the  variable, 
s  =  (R  +  H)/RN/t* 


(19) 


and  substitute  it  into  Eq.  (18)  to  give. 


After  performing  the  integration,"^  the  solution  is: 


(20) 


C  = 


(21) 


Initially,  the  concentration  within  the  chimney  was  unity.  Therefore,  the  number  of 
chimney  volumes  of  gas  expended  in  purging  is: 


N  =  1/C  . 


(22) 


Four  numerical  examples  for  H/R  values  of  2,  4,  5,  and  10,  which  shotild  cover  the 
range  of  interest  for  most  nuclear  chimneys,  are  given  in  Tables  I  and  II,  and  are 
shown  graphically  in  Fig.  1.  Other  examples  for  H/R  values  of  2.5,  2.875,  3,8,  and 


Table  I.  Numerical  computations  for  Eq.  (17b): 


C(t*)  = 

Mr) 

t* 

H/R  =  2 

H/R  =  4 

H/R  =  5 

H/R  =  10 

1 

1.000 

1.000 

1.000 

1.000 

4 

0.369 

0.569 

0.615 

0.710 

9 

0.000 

0.317 

0,387 

0.542 

16 

- 

0.102 

0.422 

25 

- 

- 

0.000 

0.328 

36 

- 

- 

0.252 

49 

- 

- 

- 

0.188 

64 

- 

- 

- 

0.133 

81 

- 

- 

- 

0.084 

100 

- 

- 

- 

0.040 

121 

- 

- 

0.000 
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Table  IL  Numerical  computations  for 
Eqs.  (21)  and  (22): 


and  N  =  1/C 


H/R 

C 

N 

0.2 

0,404 

2.48 

0,4 

2.98 

3,36 

0,5 

0.271 

3.69 

10 

0.207 

4.83 

5,0,  which  correspond  to  various  test 
assemblies,  are  plotted  in  Fig.  2  and 
Figs.  5  through  14  for  comparison  with 
both  the  uniform  mixing  model  (described 
below)  and  the  various  test  data. 

The  classic  uniform  mixing  model  is 
described  by  the  equation: 

=  exp  -  (V^/Vq)  (23) 


- where  is  the  concentration  at  any  time, 

t,  the  total  production  up  to  time,  t,,  the  chimney  volume,  and  with  an  initial 
concentration  set  at  unity.  The  decrease  in  concentration  as  a  function  of  chimney 
volume  produced  is  shown  in  Figs.  2a  and  2b,  Compared  with  concentrations  deduced 
from  the  displacement  model  (e.g.,  Eq.  17b).  The  uniform  mixing  concentrations  are 
independent  of  chimney  shape. 

These  two  models  are  compared  qualitatively  in  Figs.  3a,  b,  and  c.  Initially, 
both  are  the  same  as  shown  in  Fig.  3a,  however,  as  production  proceeds,  gas  from 
the  nearer  chimney  walls  sweeps  all  contaminated  gas  out  in  front  of  it,  successively 
purging  deeper  and  deeper  into  the  chimney.  However,  if  mixing  is  rapid  all  of  the 
gas  is  contaminated  equally  at  any  given  time,  but  the  concentration  decreases 
exponentially. 


Fig.  1.  Estimateci  radioactivity  concentration  (initial  value  =  100  percent)  during  well 
production. 
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Fig.  2a.  Comparison  of  theoretical  values  of  different  R's  to  perfect  mixing. 

In  the  above,  the  time  unit  t*  has  been  related  to  the  time  required  for  transit 
through  the  chimney  along  the  z  =  0  streamline.  The  time  required  for  removal  of 
one  chimney  volume  of  gas  has  not  been  related  except  in  an  indirect  manner  through 
C  and  the  N  that  is  obtained  by  the  integration  of  C(t*)  with  respect  to  t*  between  the 
limits  of  t*(0)  and  t*(H),  the  period  of  purging.  Now  t*(H),  Cand  N  per  Eqs.  (8),  (21) 
and  (22)  are  functions  of  R  and  H.  If  all  parameters  are  held  constant  except  H,  it 
follows  that  the  time  for  purging  is  proportional  to  H  and  N: 


t*(H)  =  pHN  (24) 

where  |3  is  a  constant  of  proportionality.  Then  from  Eq.  (24),  for  chimneys  a  and  b, 
t*(H  )  H  N 

<25) 

Then  if  t  (H),  H  and  N  are  known  for  one  case,  Eq.  (25)  can  be  used  as  well  as  Eq.  (8) 
to  find  t*(H)  for  another.  This  is  demonstrated  in  Table  111  with  H  =  5R  as  the  reference 
case. 
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Therefore,  the  mathematical  model  does  not  scale  from  one  chimney  to  another  without 
some  error,  but  this  is  not  surprising  in  view  of  the  many  approximations.  The 
chimneys  with  smallest  H/R  should  be  most  nearly  approximated  by  this  treatment. 


EXPERIMENTAL  WORK 


In  an  attempt  to  model  the  purging  of  radioactivity  from  a  nuclear  chimney,  a 
series  of  tests  were  performed  and  are  compared  with  the  theoretical  models. 
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Description 

The  apparatus  used  to  test  the  models  is  shown  in  Figs.  4  and  5.  It  consists  of  a 
set  of  concentric  cylinders  set  in  a  glass  cylindrical- walled  bowl.  The  innermost,  which 
has  a  drain-line  beneath,  represents  the  chimney  and  production  well.  For  some  test 
runs  this  inner  cylinder  was  filled  with  glass  beads;  for  other  test  runs  the  beads  were 
removed  to  permit  mechanical  stirring  of  the  solutions.  The  inner  cylinder  is  separated 
from  the  next  sand-filled  region  by  a  stainless  steel  screen  and,  during  setup,  by  a  lucite 
barrier  tube.  The  sand-filled  region  representing  the  gas  reservoir  is  separated  from 
the  outer  container  wall  by  a  perforated  lucite  cylinder  so  there  is  a  manifold  for  uni¬ 
form  flow  surrounding  the  low  permeability  sand-filled  region. 

The  height  and  diameters  of  the  innermost  cylinder  and  sand  were  different  in 
the  runs  to  simulate  various  h  to  r  ratios  of  nuclear  chimneys. 

Water  was  introduced  into  the  outermost  region  flooding  the  annulus  and  low 
permeability  sand  to  a  level  just  below  the  top  of  the  sand.  Dilute  sulphuric  acid  was 
introduced  in  the  center  cylinder  to  the  same  level.  The  central  cylinder  drain -line 
representing  the  production  well  led  through  a  flow-control  valve  and  recording  pH 
meter.  The  dilute  acid  represented  the  contaminated  gas  and  the  water  in  the  sand 
represented  the  uncontaminated  gas. 

Procedure 

Prior  to  the  start  of  each  flow  test,  the  outer  lucite  manifold  and  the  thin -walled 
cylinder  were  placed  in  the  bowl.  Beads  were  added  in  the  small  center  cylinder,  and 
sand  in  the  space  between  it  and  the  outer  lucite  cylinder.  Water  and  acid  were  then 
added. 


Fig.  4.  Cross-section  of  the  apparatus  used  in  testing  the  models. 


Fig.  5a.  The  above  apparatus  was  used  to  test  the  models.  Note  that  the  lucite  dia¬ 
phragm  and  center  cylinder  or  chimney  is  full  of  beads.  (No  sand  is  present.) 
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lucite  diaphragm  and  center  cylinder  with  the  stainless  steel  screen  holds  the  sand  in  plac 
r  was  removed  at  the  start  of  the  run. 


Fig,  5c.  The  experimental  set-up,  showing  the  sand-filled  bowl  at  top  with  the  pH 
meter  and  recorder,  and  the  flowmeter  and  timer  at  the  bottom. 
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The  flow  rate  through  the  micro-valve  was  adjusted  and  set,  the  pH  meter  was 
calibrated,  and  the  water  and  acid  levels  were  adjusted  at  about  10  mm  below  the 
sand-bead  level,  to  help  prevent  "short  circuits"  across  the  top  of  the  sand.  A  timer 
was  started,  the  ON-OFF  valve  opened,  and  the  thin-walled  barrier  cylinder  between 
the  "chimney"  and  production  sand  gently  removed. 

The  gradual  change  in  pH,  from  the  initial  pH  of  the  starting  acid  used  in  the 
center  cylinder  to  the  pH  of  the  water  in  the  outer  area,  was  recorded  on  a  timed  chart. 
A  total  of  14  tests  were  run  and  the  data  presented  in  Table  IV  and  Figs.  6  through  14. 


Table  IV.  Summary  of  flow  test  data. 


Liquid 

Cylinder 

Average 

height 

radius 

Ratio 

flow 

Volume 

_ pH _ 

R\in 

h 

r 

'  R 

rate 

Acid 

Water 

Acid 

Water 

Notes 

No. 

(mm) 

(mm) 

(h/r) 

cc/min 

cc 

cc 

t*^  (water-acid)22‘^C 

1 

100 

40 

2.5 

>600.0 

150^ 

NA 

0.2 

NA 

Coarse  sand. 

2 

100 

40 

2.5 

-100.0 

150 

1100 

0.2 

8.7 

20  mesh,  no  mani¬ 
fold,  fast  flow 

3 

100 

40 

2.5 

-'lOO.O 

160 

1200 

2.0 

6.5*^ 

.  rates. 

4 

115 

40 

2.875 

170.0 

180 

1775 

1.7 

8.4 

[  Fine  sand,  -100 

5 

115 

40 

2.875 

92.0 

180 

1770 

1.7 

8.4  j 

mesh,  with  Incite 

1  manifold. 

6 

100 

40 

2.5 

42.0 

170 

1550 

1.7 

8.4 

7 

152 

40 

3.8 

21.5 

195 

2500 

1.7 

8.4 

8 

152 

40 

3.8 

17.0 

195 

2500 

1.66 

„  f  Start-up  trouble 

(seal) 

9 

192 

40 

4.8 

20.0 

270 

3000 

1.82 

8.9  j 

[  Aborted,  due  to 

1  leak. 

10 

175 

35 

5.0 

24.0 

260 

3700 

2.25 

8.9 

11 

175 

35 

5.0 

16.0 

260 

3750 

2.00 

8.9 

12 

175 

35 

5.0 

17.0 

265 

3750 

1.42 

8.9 

1  No  beads  in  center 

13 

160^^ 

32 

5.0 

26.0 

470 

2780 

1.58 

6.2^  1 

,  column  A  larger 
column  vol. 

1 

1  Run  14  with  stirrer 

14 

160 

32 

5.0 

26.0 

450 

2750® 

1.55 

6.5  1 

in  center  column, 

’  Constant  mixing. 

Plotted:  4,  5,  8,  12,  13,  14  and  theoretical 


a 

Drainage  vol.  Total  acid  on  to  dry  beads  175. 

After  sitting  12  min.  Sand  needed  more  washing.  Later  runs  showed  little  pH 
change,  even  after  the  system  sat  for  hours. 

Q 

Well  below  sand  level  of  190  mm. 

^De -ionized  water. 

0 

Of  the  4550  cc  originally  added  to  dry  sand,  1800  cc  was  held  up  in  the  wet  sand 
after  draining  and  standing  for  10  min. 


Column  volumes 

Fig.  6.  Results  of  the  R  =  2.875  flow  test. 


Fig.  7.  Results  of  the  R  =2.5  flow  test. 


Column  volumes  Column  volumes 

Fig.  8.  Results  of  the  R  =  3.8  flow  test.  Fig.  9.  Results  of  the  R  =  3.8  flow  test 

(slow  rate)  based  on  the  Higgins - 
Rodean  fluid-displacement  theory. 


Fig.  10.  Results  of  the  R  =  5.0  medium 
fast -flow  test. 


Fig.  11.  Results  of  the  R  =  5.0  slow-flow 
low -acid  flow  test. 


Fig.  12.  Results  of  the  R  =  5.0  very  slow- 
flow  test  with  high-acid. 


Fig.  13.  Results  of  the  R  =  5.0  very  slow- 
flow  test.  (No  beads,  larger 
chimney  volume,  purer  water, 
and  high  acid.) 
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The  acid  concentration  determined 
from  pH  readings  represents  the  radio¬ 
activity  concentration  and  its  decrease 
simulates  the  decrease  in  radioactivity 
during  gas  production. 

Discussion 

The  first  model  tested  in  runs  1 
through  12  is  "simple  fluid -displacement" 
without  diffusion  or  dispersion.  With  this 
model  there  is  no  perceptible  change  in  the 
concentration  in  fluid  first  leaving  the 
chimney.  Then,  as  the  flow  continues,  the 
concentration  decreases  fairly  rapidly,  as 
indicated  in  Figs.  6  through  12,  which 
depict  the  theoretical  model  for  displace¬ 
ment  and  relate  C/Cq  to  column  volume. 

In  a  finite  time,  zero  concentration  is 
reached  and  the  chimney  is  completely 
purged.  In  the  simplified  theory,  it  is 
assumed  that  there  is  a  piston-like  dis¬ 
placement  with  no  dispersion  or  diffusion 
along  each  infinitesimal  line.  However,  in 
view  of  the  complex  geometry  followed  by  the  stream  lines  from  the  walls  of  the  cylinder 
to  a  point  in  one  end  of  the  cylinder,  the  concentration  varies  with  time  in  a  complex 
way.  [See  Eqs.  (17a)  and  (17b).] 

These  experiments  show  that  the  concentration  in  the  discharge  remains  relatively 
high  at  first,  and  then  begins  to  drop  fairly  rapidly.  This  is  qualitative,  if  not 
quantitative  agreement  with  the  simple  fluid -displacement  model  theory.  However,  the 
concentration  is  significantly  above  zero  at  late  times.  The  chimney  has  not  been 
completely  purged  and  the  model  tends  to  resemble  the  perfect  mix  situation,  in  place 
of  the  simple  fluid -displacement  model.  Mixing  becomes  evident  in  the  later  stages 
because  of  the  dispersion  and  diffusion  that  has  occurred  and  because  the  theory  is  only 
approximate.’ 

In  Figs.  6  through  14,  the  change  in  acid  concentration  versus  the  cumulative  flow 
measured  in  chimney  (inner  cylinder)  volumes  for  each  of  the  test  runs  and  is  to  be 
compared  with  the  solid  lines  from  the  calculations  best  fitting  and  test  situation.  The 
dashed  lines  in  Figs.  13  and  14  are  for  perfect  mixing  in  the  chimney. 

Run  4  was  an  anomaly,  because  the  integral  under  the  curve  was  much  greater 
than  one  (one  full  chimney  volume,  over  C/Cq  range).  (Evidently,  an  error  in  pH  meter 
calibration  occurred  during  this  run.)  Runs  1,  2,  and  3  were  used  to  adjust  the 


Theoretical  fluid 
displacement  theory  ^ 
R  =5.0 


Theoretical 
perfect  mix 


1  2  3  4  5 

Column  volumes 


Fig.  14.  Results  of  the  R  =  5.0  very  slow- 
flow  test  "perfect  mix."  (No 
beads,  larger  chimney  volume, 
and  constant  stirring.) 


-18- 


mechanics  of  the  system  and  run  9  was  aborted  because  of  a  leak.  R\ins  7  and  8  are 
very  similar  in  R  and  flow  rate,  as  are  runs  11  and  12.  The  differences  between 
members  of  these  pairs  indicate  the  reproductibility  of  the  experimental  system. 

In  the  second  theoretical  model,  tested  in  run  No.  14,  there  is  continuous 
"perfect  mixing"  during  production  flow.  Under  these  conditions,  the  fluid  leaving  the 
chimney  immediately  begins  to  show  the  effect  of  the  fluid  flowing  into  the  chimney 
from  outside.  Theoretically,  the  chimney  will  never  be  purged  completely  because 
there  is  always  some  of  the  original  chimney  material  present;  the  concentration  of 
original  material  gradually  approaches  zero,  but  is  equal  to  zero  only  at  infinite  time. 
The  results  shown  in  Fig.  14  are  very  close  to  the  theoretical  model. 

Run  13  (see  Fig.  13)  was  performed  without  glass  beads  and  without  stirring  in 
the  center  volume  in  order  to  determine  the  maximum  amount  of  diffusion  and  mixing 
which  was  caused  by  the  acid  concentration  gradient.  It  can  be  observed  that  the  re¬ 
sults  resemble  rxons  11  and  12  closely,  so  the  mixing  induced  by  acid  gradients  (while 
obviously  present)  probably  is  not  perturbing  results  grossly. 

In  summary,  the  results  indicate  the  experimental  data  falls  between  the  two 
theoretical  extremes:  (1)  conditions  resemble  simple  fluid -displacement  flow  at  the 
start  and,  (2)  at  late  times,  diffusion  and  dispersion  result  in  a  condition  more  nearly 
related  to  "perfect  mixing."  The  test  in  which  stirring  and  constant  mixing  were 
employed  (run  14),  closely  corresponds  to  the  theoretical  curve  and  is  probably 
indicative  of  the  general  accuracy  of  the  test  method. 

Conclusions 

It  can  be  concluded:  (1)  in  the  absence  of  mechanical  mixing,  the  fluid  in  the  more 
permeable  central  column  was  purged  more  nearly  in  accordance  with  the  "displacement" 
theory,  (2)  the  experimental  setup  was  not  perfect,  in  that  flow  through  the  fine  sand 
could  be  variable  over  the  surface  of  the  cylinder,  (3)  it  is  likely  that  production  of 
radioactive  gas  from  nuclear  chimneys  in  gas -bearing  rock  will  be  influenced  by 
several  processes  such  as:  (a)  convective  mixing  leading  to  more  nearly  exponential 
decay  of  effluent  concentrations;  (b)  nonuniform  production  from  the  chimney  walls; 
and  (c)  intermittent  rates  of  production  required  by  other  constraints.  All  of  these 
factors  will  tend  to  make  the  real  field  observations  fall  between  the  extremes  presented 
by  the  two  theoretical  models. 

The  tests  are  sufficiently  convincing  to  warrant  purge  tests  of  nuclear  chimneys 
when  it  becomes  desirable  that  all  the  radioactive  gases  be  reduced  100 -fold  or  more  in 
concentration  with  minimum  gas  flaring. 

In  order  to  minimize  this  total  amount  of  gas  flared  during  purging,  the  purging 
rate  (flow -rate)  should  be  as  high  as  possible  to  minimize  thermal  convective  mixing 
and  the  purging  should  proceed  uninterruptedly.  Any  interruption  will  allow  mixing 
to  recontaminate  those  areas  already  cleaned  by  incoming  gases.  It  is  unfortunate  that 


field  tests  such  as  those  conducted  on  Project  Gasbuggy  prior  to  the  fall  of  1968  were 

85 

interrupted  and  provided  almost  no  new  data.  From  the  Kr°^  concentration  reported 

5 

by  Smith  and  Momyer  it  can  be  concluded  that  initially  the  concentrations  were  con¬ 
stant,  consistent  with  the  displacement  model,  but  that  convective  mixing  was  occur¬ 
ring  at  such  a  degree  that  after  shut-in  of  several  days  duration,  the  contaminant 
would  be  uniformly  remixed. 
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